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ABSTRACT 

We present the 2-60 keV spectrum of the supernova remnant Gassiopeia A 
measured using the Proportional Gounter Array and the High Energy X-ray 
Timing Experiment on the Rossi X-ray Timing Explorer satellite. In addition 
to the previously reported strong emission-line features produced by thermal 
plasmas, the broad-band spectrum has a high-energy "tail" that extends to 
energies at least as high as 120 keV. This tail may be described by a broken 
power law that has photon indices of Fi = LSt^g and F2 = 3.041qj| and a 
break energy of = l5.9toi keV. We argue that the high-energy component, 
which dominates the spectrum above about 10 keV, is produced by synchrotron 
radiation from electrons that have energies up to at least 40 TeV. This conclusion 
supports the hypothesis that Galactic cosmic rays are accelerated predominantly 
in supernova remnants. 

Subject headings: ISM: individual (Gassiopeia A) — supernova remnants — 
radiation mechanisms: non-thermal — acceleration of particles — cosmic rays 
— X-rays: general 
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1. Introduction 



Galactic cosmic rays that have energies below ~ 10^^ eV are believed to be accelerated 
predominantly in the shocks of supernova remnants (SNRs). Prior to the last two years, 
little experimental evidence existed to support this hypothesis. Radio measurements 
show that many SNRs emit synchrotron radiation from nonthermal electrons, but these 
measurements are not sensitive to emission from electrons that have energies greater than 
~ 10^'^ eV. Gamma-ray measurements at energies < 10^° eV also show that some SNRs 
contain nonthermal particles ( [Esposito et al. 1996|) , but the interpretation of these results is 
controversial ( [Allen et al. 1995| ; [Lessard et al. 1995| ; |de Jager fc Mastichiadis 1997| ; [Gaisser 



Protheroe, fc Stanev 199^ ; [Sturner et al. 19'97 ). However, recent measurements of X-ray 
spectra of SNRs have provided the best evidence to date that SNRs may be responsible for 
most of the cosmic-ray acceleration in the Galaxy. 

While the X-ray spectra of SNRs have traditionally been modeled as thermal 
emission from shock-heated plasmas, the spectra of SN 1006 QKoyama et al. 1995|) and 
RX J1713.7— 3946 ([Koyama et al. 1997| ) are dominated by nonthermal emission interpreted 
to be synchrotron radiation from electrons having energies up to ~ lO^'^ eV. IC 443 may 
also contain a site that emits X-ray synchrotron radiation from electrons having energies at 
least as high as 2 x 10^^ eV ( [Keohane et al. 19ir7| ). 



We present the 2-60 keV spectrum of Cassiopeia A (Gas A) measured using detectors 
on the Rossi X-ray Timing Explorer (RXTE) satellite. This spectrum has a high-energy 
"tail" that extends to energies at least as high as 60 keV. The strong emission lines in the 
0.5-10 keV spectrum ( [Holt et al. 1994| ; [Jansen et al. 198^ ; [Isunemi et al. 1986| ) can be 
described by two-temperature thermal plasma models whose flux is dominated by emission 
from either swept-up material ( [Borkowski et al. 1996| ) or shocked ejecta ( |Vink, Kaastra, S 
Bleeker 1996| ; [Jansen et al. 198"BD . However, these models do not describe the higher-energy 
tail of the spectrum (l^'avata et al. 1997; Pravdo & Smith 1979[ ; [Hatsukade & Tsunemi l9^ ) 
that has been detected up to at least 120 keV ( P?he et al. 1996[) . The results of our analysis 
imply that both the X-ray tail and the radio spectrum {Su = 2723 (z//l GHz)~° '''^ Jy, 
Paars et al. 1977|) are produced by synchrotron radiation from a common population of 
nonthermal electrons whose spectrum extends up to at least 4 x 10^^ eV. 



2. Data and Analysis 



Between 1996 March 31 and 1996 April 17, Gas A was observed for 186 ks using 
the Proportional Gounter Array (PGA) and the High Energy X-ray Timing Experiment 
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(HEXTE) on the RXTE satellite. The PC A ( |Jahoda et al. 1996|) and HEXTE ( |Gruber et aT 



|1996| ) comprise a 2-250 keV spectrophotometer that has an energy resolution ~ 10-30% 
and a collecting area of about 7000 cm^ at 6 keV and of 800 cm^ at 60 keV. 

Figures I and | show the PCA and HEXTE spectra of Cas A for 91 ks of the PCA 
data and for 96 ks of HEXTE data. Because the strong emission lines may be fit by 
models with two thermal plasmas and with fluorescent Ka emission from iron in dust 
grains ( [Borkowski fc Szymkowiak 1997| ), we model these features using two Raymond-Smith 
components and a 6.40 keV Gaussian. While such a model provides an adequate fit to the 
0.5-10 keV spectrum, it does not fit the spectrum above ~ 10 keV (fig. At least one 
more component is needed to model the broad-band X-ray spectrum. 

We fit the broad-band spectrum using models that include only one additional 
component that is either a thermal bremsstrahlung, a power law, a power law that has an 
exponential cutoff, or a broken power law. In general, the high-energy spectrum is better fit 
by the nonthermal components than a thermal bremsstrahlung component. Furthermore, 
the spectrum is better fit by a broken power law than a power law or a cut off power law. 
Table |1| summarizes the best fit models that include either a broken power law component 
or a thermal bremsstrahlung component. Abundances of some line-emitting elements are 
listed in the table because the abundances depend on the high-energy component. The 
abundances should be regarded as rough estimates because a better treatment requires 
plasmas that are not in ionization equilibrium. The values of in table || were computed 
for only the > 10 keV portion of the PCA and HEXTE spectra where the high-energy 
component dominates the flux. The "la" errors quoted in this paper, which may be 
underestimated, were computed by finding the extrema of the contour that has a value 
of which is larger than the minimum value by 2.3. 

Figure |I| shows the spectrum of the best fit model "folded" through the PCA and 
HEXTE response matrices. The lower panel of this figure shows the ratio of the measured 
spectra to this model. Aside from the uncertainty in the cross calibration of the PCA 
and HEXTE, this ratio is dominated by systematic errors in the response matrix and the 
background spectrum for the PCA and by statistical errors for the HEXTE. For example, 
the features near 2 keV and between 5 and 8 keV in the lower panel of figure reflect 
inaccuracies of the v2.0.2 PCA response matrix that are observed when modeling PCA 
spectra for other sources, such as the Crab ( jJahoda et al. 195^ . These features are not 



evident in the fit to the spectrum of the Gas Imaging Spectrometer no. 2 (GIS2) on the 
Advanced Satellite for Cosmology and Astrophysics. 

Figure | shows the GIS2, PCA, HEXTE, and OSSE (Oriented Scintillation 



Spectrometer Experiment on the Compton Gamma Ray Observatory satellite, [The et al 
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1996) spectral data and the "unfolded" spectrum of the best fit model (solid curve). 
For comparison, figure ^ also shows the spectrum of a model that includes only the two 
Raymond-Smith components and the Fe Ka dust component (dashed curve). 



3. Discussion 

The spectra of several emission mechanisms — thermal and nonthermal bremsstrahlung, 
a pulsar, inverse Compton scattering, and synchrotron radiation — could be described by a 
broken power law in the range of 10-63 keV. Each possibility is discussed in the context of 
Cas A. 

It is unlikely that the high-energy X-ray tail of Cas A is thermal. As shown in 
figure and table |l|, the broad-band X-ray spectrum is not well fit by two or three thermal 
component models. Furthermore, a physical interpretation of a model that has three 
thermal components is difficult. Two of the components could be associated with the 
forward-shocked circumstellar material and the reverse-shocked ejecta. The third might 
result from a nonuniform distribution of the circumstellar or ejected matter. However, 
spatially-resolved spectroscopy shows that Cas A has similar plasma conditions across the 
remnant ([Holt et al. 1994|) . 



Asvarov et al. (1989) suggest that the high-energy X-ray spectrum of Cas A 
is dominated by nonthermal bremsstrahlung emission. However, an estimate of the 
nonthermal bremsstrahlung spectrum (fig. ^ yields a 20-50 keV effective photon index of 
~ 1.9, which is significantly smaller than the fitted index of 3.04loj|. Therefore, the X-ray 
tail of Cas A appears to be too steep to be consistent with a nonthermal bremsstrahlung 
emission process. 

Although pulsars can produce nonthermal X-ray continua, it is unlikely that the 
high-energy tail of Cas A is produced by a pulsar. Analyses of X-ray ( |Allen et al. 1997| ; 



Fravdo fc Smith 1979| ), radio ( |Woan fc Duffett-Smith 1993| ) and optical ( [Horowit"z^ 



Fapaliolios, fc Carleton 1971| ) data reveal no evidence of pulsation and no bright point-like 



feature is observed near the center of the remnant ( [b'abian et al. 19ED| ; [Holt et al. 1994 
Jansen et al. 1988|) . 



An inverse Compton spectrum could be consistent with the shape of the tail, but the 
estimated flux is a factor of ~ 10^ smaller than the measured flux at energies in the range 
of 10-63 keV (flg. |). 



In contrast to the other emission processes, both the shape and the flux of the 
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high-energy component of the X-ray spectrum of Cas A are consistent with a synchrotron 
radiation mechanism. In general, the synchrotron spectrum of a SNR is expected to 
span the entire range from radio to X-ray energies and to steepen gradually at or near 
X-ray energies ( Reynolds 1997 ). This expectation is consistent with the multiwavelength 



spectrum of Cas A because the extrapolation of the radio synchrotron spectrum to X-ray 
energies does not lie above the observed thermal emission at infrared wavelengths or below 
the X-ray spectrum (fig. |^) and because the broken power law of table |I] approximately 
describes a gradually steepening spectrum. The synchrotron spectrum of figure ^ is a 
power-law spectrum with a e~^^'^^^ ^'^^ cutoff. This photon spectrum corresponds to a 
power-law electron spectrum with a Q-i^<=/^'>^'^^^'^ cutoff and assumes that the spectral 
distribution of synchrotron photons emitted by electrons of a common energy is a delta 
function. A better estimate of the synchrotron spectrum, which is expected to steepen more 
gradually than g-V-^^/i ^ev ([Reynolds 1997| ), requires a numerical simulation of the particle 
acceleration conditions appropriate for Cas A including the effects (1) of synchrotron losses 
on the electron spectrum, (2) of possible curvature in the electron spectrum ( [Mezger et alj 
1986| ; [Ellison fc Reynolds 1991| ), and (3) of the spectral distribution of synchrotron photons 



for electrons of a given energy. 



Figure ^ shows the radio ( Paars et al. 197"7|) , infrared ([Mezger et al. 1986|) , RXTE 



OSSE ([I'he et al. 19961) , EGRET ( [Esposito et al. 199q ), and Whipple QLessard et al. 



1995|) spectral data of Cas A. The estimated synchrotron spectrum is a single power law 



that has a e~ v^^^^^ cutoff. Also shown are estimates of the photon spectra produced 
by nonthermal bremsstrahlung, by inverse Compton scattering of the cosmic microwave 
background radiation, and by the decay of neutral pions. These latter three estimates were 
computed for electron and proton spectra that have a common spectral index of 2.54 using 
equation 13 and figure 3 of Gaisser et al. (1997). The electron density spectrum is assumed 
to be 7 X 10~^ cm~^ GeV~^ at 1 GeV. The density of accelerated nuclei is assumed to be 
one hundred times larger than the density of accelerated electrons independent of energy. 
The nonthermal bremsstrahlung and 7r° decay spectra were computed for a density of nuclei 
in Cas A of 30 cm~^. The shape of the nonthermal bremsstrahlung spectrum at energies 
below 1 MeV is based on the shape of the nonthermal bremsstrahlung spectrum of Sturner 
et al. (1997, fig. 11). 



4. Conclusion 



The RXTE data reveal that the X-ray continuum of Cas A has a nonthermal high-energy 
"tail." The tail (1) is qualitatively consistent with a simple model of synchrotron emission 
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from SNRs, (2) is inconsistent with the expected shape of a nonthermal bremsstrahlung 
spectrum, and (3) is inconsistent with the estimated inverse Compton flux. No evidence of 
pulsation or of a bright central point-like feature has been reported. Therefore, the tail is 
most likely produced by synchrotron radiation. 

X-ray synchrotron emission has important implications for particle acceleration in 
Cas A. The energy of a synchrotron photon, E^, is related to the energy of the emitting 
electron and the magnetic field strength by i?^ ~ 0.6i?^G-^i4 keV, where B^q is the 
magnetic field strength in units of /xG and £^14 is the energy of the electron in units of 
10^^ eV. If the synchrotron spectrum extends to energies at least as high as 120 keV and 
the magnetic field is 1 mG (the "equipartition" value, e.g. [Longair 1994| ), the accelerated 



electron spectrum extends up to at least 4 x 10^^ eV. Since electrons and protons are 
expected to be accelerated in the same manner at energies much larger than the rest mass 
energy of the proton ( Ellison fc Reynolds 1991 ), this result provides strong evidence for 



the acceleration of cosmic-ray protons in Cas A to the same energies. The calculation of 
the equipartition value of the magnetic field leads to an estimate of the total amount of 
energy in accelerated protons and electrons of ~ 3 x 10^^ erg ( [Longair 1994 ). This amount 



of energy is comparable to the average amount of energy required from a Galactic SNR 
(~ 3-10 X 10'^^ erg), over the lifetime of the remnant, if all of the Galactic cosmic rays are 
accelerated in SNRs. 

Cas A is the fourth SNR reported to exhibit evidence of X-ray synchrotron radiation. 
Collectively, the resuhs for Cas A, SN 1006 (|Koyama et al. 1995|) , RX J1713. 7-3946 
( [Koyama et al. 1997|) , and IC 443 ( [Keohane et al. 1997| ) support the hypothesis that 



Galactic cosmic rays are accelerated predominantly in SNRs. 
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Fig. 1. — The PCA and HEXTE count spectra of Cas A. The histogram is the broken 
power-law model of table |l[ The lower panel shows the ratio of the data to the histogram. 
See text for details. 



Fig. 2.— The GIS2, PCA, HEXTE, and OSSE (The et al. 1996) photon spectra of Cas A. 
The solid curve shows the broken power-law model of table 0. The dashed curve shows what 
the shape of the spectrum is if no high-energy component is included in the fit. 



Fig. 3. — The multiwavelength photon spectrum of Cas A. The four broken curves are 
estimates of the fluxes from synchrotron radiation (S), from nonthermal bremsstrahlung 
(NB), from inverse Compton scattering of the cosmic microwave background (IC), and from 
the decay of neutral pions (vr°). See text for details. 
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Table 1. Parameters of the Spectral Models'^ 



Broken Power Law 










Abundance'' 






{E < E^) 


Eb 
(keV) 


r2 

{E > Eb) 


Flux at 1 keV Si 
(ph cm-2 s-i keV-i) 


S Ar Ca 
(Solar) 


Fe 


(69 dof) 




-"-■^■■'-(1.4 




0.038 2.1 


4.3 8.3 2.0 


0.0 


105 


Thermal Bremsstrahlung 




kT 






Abundance'' 








(keV) 




(10^* cm-3) 


(as above) 




(71 dof) 








1.16 ±0.04 2.0 


4.1 7.7 2.3 


0.6 


236 



'^These fits also include two Raymond-Smith components that have temperatures of kTi = 0.6 keV 
and kT2 = 2.9 keV and a 6.40 keV Gaussian component (1.3 x 10~^ photons cm~^ s~') for Ka 
emission from iron in dust grains. The fit column density rin = 1-1 x 10^^ H atoms cm^^. 

''The Si, S, and Ar abundances pertain to the kTi = 0.6 keV plasma. The Ca and Fe abundances 
pertain to the kT2 = 2.9 keV plasma. 

"^For the energy range from 10 to 63 keV. 

''The quantities Ue and rii are the electron and ion densities respectively. The average value of 
neTii is ~ 25 cm~^ for a distance of 3.4 kpc and a filling fraction of 1/4. 
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